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Using a general Green function formulation, we re-derive, both, (i) Spitzer and his followers results 
for the winding angle distribution of the planar Brownian motion, and (ii) Edwards-Prager-Frisch 
results on the statistical mechanics of a ring polymer entangled with a straight bar. In the statistical 
mechanics part, we consider both cases of quenched and annealed topology. Among new results, we 
compute exactly the (expectation value of) the surface area of the locus of points such that each of 
them has linking number n with a given closed random walk trajectory (= ring polymer). We also 
consider the generalizations of the problem for the finite diameter (disc-like) obstacle and winding 
within a cavity. 
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I. INTRODUCTION 

In 1958, Spitzer [jj established the following result. 
Consider the two-dimensional random walk starting at a 
point other than 0, and let 9{t) be the total continuous 
angle wound by the walker around O up to time t (see fig- 
ure^a). The Spitzer law says that the quantity 9{t)/\ut 
at large enough t is Lorenz (or Cauchy) distributed: 
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With a remarkable delay of about 25 years, a large group 
of followers studied this law in depth 11 H lEflaB ll 
S IIS IH 13 111 The central finding of these studies 
attributes the divergent moments of the Spitzer distribu- 
tion Q, e.g. (0 2 ), to the small scale properties of the 
regular random walk trajectories. Simply speaking, in- 
finitely large winding is accumulated while the trajectory 
is wandering infinitely close to the obstacle O. Accord- 
ingly, this pathology of divergent moments is removed by 
incorporating any kind of "granularity," or short length 
scale cut-off, in the model. Such modification of the 
model can be achieved in quite a few ways. One way 
is to consider the random walk on the lattice instead 
of the continuous space ;9j; another way is to look at 
the winding around a finite obstacle, say, a disc of some 
radius b UE1 0! one more possibility is to examine a 
broken line of straight segments of finite length b each 
instead of standard Wiener-measured random walk; yet 
another way is to consider a worm-like smooth curve with 
an effective segment b (that is, the curve which adopts 
smoothly curved shapes r(s) with the weight propor- 
tional to exp [—(6/2) f r 2 ds] , where s is the arc length). 
In all of these cases, winding is characterized by the non- 
pathological distribution 
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A similar distribution is also characteristic for the wind- 
ing of the self-avoiding walk UlEHi self-avoidance, in 
this case, is just another way to suppress infinite wind- 
ing at infinitely small length scale. Mathematically, it 
turns out that the winding angle distribution is in fact 
an example of a broad class of limiting laws for the two- 
dimensional random walk 0, H . 

Many studies of winding angle distribution 0, 0, H IE 
0, E3, [3 El claim that entanglement of long polymer 
filaments is (one of) their motivation(s). Indeed, the re- 
lation to polymer physics does exist. It was found in 
1967, almost a decade after Spitzer by Edwards Q 
and, independently, by Prager and Frisch [H^ (see also 
an influential review [l6j)- These authors came up with 
the model of a polymer chain wound around a straight 
bar and topologically entangled with this bar. Given the 
analogy of a polymer chain conformation with the ran- 
dom walk trajectory, the Edwards-Prager-Frisch model 
is essentially the same as that examined by Spitzer p|. 
Neither of the works 0, fl5l H(| makes a reference to 
|j| . Most likely, mathematical work was not known to 
physicists at the time, but even apart from that, authors 
of the works UJ,h5lll6| did not examine winding angle 
distribution for the random walk with open ends, their 
goal was obviously to compute quantities similar to those 
of physical interest for real polymers. Unfortuna tely , no 
explicit formula was obtained in the works 0, \loL Il6| 
comparable in simplicity to Eq. (Q. 

To our surprise, we found that this fairly old area lacks 
both unityand clarity. The studies of winding angle dis- 
tribution H, H H H|l [IHH [H contain no hint on the 
lessons of this exactly solvable model to polymer physics. 
Drossel and Kardar 0| as well as Samokhin |l2l Il3| 
brought the subject to a new level of complexity, they ex- 
amined winding angle distribution for the random walks 
in a disordered medium. Drossel and Kardar also 
provided simple derivation of the results (11121) and applied 
it to many physical situations involving directed poly- 
mers, but all that yields little insight into the topological 
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who took advantage of the more 
recent mathematical achievements 0, & IE IE IE IE E3 ■ 
Meanwhile, an exactly solvable model in general is use- 
ful if only it yields some insight (s). Upon a closer look 
and re-examination of the literature, we found that the 
model of winding can be made to meet this criteria, but 
it has not been done yet. Our plan in this paper is to re- 
consider the problem from a single common view point, 
including both winding angle distribution and some more 
physical aspects. 

Our additional motivation arises from the fact that the 
study of topological constraints in polymers in the years 
and decades after the works Hi [S Eg had been do mi- 
nated by the phenomenological approaches based on the 
reptation theory 0, At the same time, a break- 

through in microscopic understanding of this subject has 
not been achieved, and, therefore, the need for exactly 
solvable models remains high. Moreover, apart from net- 
works, there is now another large "consumer" for polymer 
topology, this is DNA physics. The DNA double helix is 
frequently found in a closed loop form, it forms knots 
of various kinds 0, |5(j , and there are special enzymes 
spending energy to simplify the entanglements pH . 

One of the key aspects of polymer topology is that 
there are two types of questions one can ask, correspond- 
ing to annealed and quenched topological disorder, re- 
spectively . The beauty of the winding model, which 
so far seems to remain underappreciated, is that it allows 
both types of questions: 

• The typical annealed topology question is that 
about ring closure experiment and knot probabil- 
ities [2E HE I27L HE EE] : having a linear polymer 
with "sticky" ends, what is the probability to ob- 
tain a certain type of a knot upon first meeting of 
the two ends [l3. 120I]? A similar question for the 
winding model is this: what is the probability that 
a random walk on the plane links number n (or 
winding angle 2im) with an obstacle? 



• The typical quenched topology question is about, 
e.g., the size or other properties of a polymer hav- 
ing a given fixed topology (e.g, knot type) [23ll2^|: 
this is necessary, e.g., to understand the diffusion 
of knotted DNA in solution or in a gel. A s imi- 
lar question for the winding model is this: given a 
polymer with fixed linking number n, what is the 
(root-mean-squared) average distance of an arbi- 
trary point on the trajectory from 0? 

To conclude the introduction, we should also men- 
tion that the shortcomings of the Edwards-Prager-Frisch 
model are well understood j^l). Basically, this model 
assumes that entanglements algebraically commute with 
each other, while the real physical situation is non- 
Abelian. 

This paper is organized as follows. In section IIIII 
we discuss the Green function formulation of the prob- 
lem and derive basic equations for all models - winding 






FIG. 1: Schematic representattion of the model, (a) Random 
walk winding around an obstacle O. This obstacle might be 
just a point, or it might be a disc of a finite radius b. (b) 
Closed polymer winding around an obstacle. Mathematically, 
this is similar to (a), except both ends are kept together, (c) 
Similar to (a), except the trajectory cannot leave a "cavity" 
of some radius B. 



around the point, around the disc, or inside the cavity. In 
section llVl we show how to re-derive and generalize the 
results (I) and ln section we consider the closed 
loop polymer, which is the random walk with connected 
ends. In section IvTl we make a few final comments. 

Our work is heavy on calculations, even though some of 
the less important ones are relegated into Appendices. As 
readers, we don't like such heavy papers. This is why we 
start from section[n]which provides an overview of major 
steps and the results for those readers not interested in 
details. 
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II. BIRD'S EYE VIEW OF THE RESULTS: FOR 
THE LAZY READER WHO DOES NOT WANT 
TO DWELL ON THE CALCULATIONS 

If you, our reader, do not want to follow our calcula- 
tions, this section offers a tour of the results for you. 

To begin with, section IIIII contains no results: it de- 
scribes the standard diffusion equation and bilinear ex- 
pansion of its Green function over the appropriate set of 
Bessel functions. Here, for a "tourist," all that is neces- 
sary to know is the notation a 2 /4 adopted for the diffu- 
sion coefficient, which means that the root-mean-square 
distance traveled by a walker during the time t is equal 
to a\ft. 

Formula (|14|l in section llV Al is our first small result, 
it is a very mild generalization of the Spitzer formula 
which takes explicit account of the distances r and r' 
of a polymer (or random walk trajectory) ends from the 
origin (or obstacle) O. Formula l|14l) gives the probability 
distribution of winding angle 9 for the random walk of 
length t with r and r' fixed. Like the Spitzer law, it has 
diverging moments, such as (9 2 ). 

The very cumbersome formula l|23ll gives a similar re- 
sult for the winding around a disc of a finite radius 6. It 
generalizes formula @ by keeping explicit track of po- 
sitions r and r' of both ends. Just like J2Jl, it decays 
exponentially and yields finite values for all moments, 
e.g., (0 2 ). One utility of this result is the analysis of 
cross-over between winding around a point with infinite 
(9 2 ) and winding around a disc with finite (0 2 ). As we 
show in Section TlV CI when the disc size 6 goes to zero, 
there opens a wide range of times t (see Eq. ffify ) where 
the probability behaves as ~ I/O 2 up to 9 about In (a/6), 
and only at larger 9 exponential decay takes over; there- 
fore, when we say that (9 2 ) diverges, this really means 
(9 2 ) ~ (In (a/6)) 2 at small 6. 

In Section |V| we go closer to the polymer view on the 
subject. For this, we consider that the two ends of the 
random walk are glued together, so that r = r' and 
9 = 27m, where (positive or negative) integer n is the 
linking number, the number of turns the polymer ring 
makes around the obstacle. Figure |21 depicts the statis- 
tical weight of the polymer conformations with the link- 
ing number n as a function of r 2 /t. Qualitatively, this 
exhibits a behavior similar to that of the knotting prob- 
ability as a function of chain length, because all cases 
with n / (similar to non-trivial knots) reach maximal 
weight at some intermediate values of r and/or t. 

Sections IV Bl and IV CI present our most original and 
most interesting findings. In particular, we consider 
the following question. Given the closed random walk 
trajectory of the length t, we consider £„ - the locus 
of points around which the trajectory makes exactly n 
turns. Then, what is the area of £„? We denote by 
a n (t) the average (over the random walk trajectories) of 
this area, then formula l|44l) provides the exact answer to 
this question. The essence is that <r n (t) decreases very 
slowly with n, only as 1/n 2 . Of course, ideologically this 



is similar to the slow decay of the Spitzer distribution 
Another look at the same result is to think about 
a virial coefficient of a polymer ring with a long straight 
bar. Their interaction is topological in nature |25j |. and 
the virial coefficient can be understood as the surface ex- 
cluded for a polymer ring by the presence of the obstacle 
if the ring is not entangled with the obstacle. This virial 
coefficient is the sum of all cr n (t) with n / , and it is 
exactly equal to 7rfa 2 /f2. 

Note that the former view of a n (t) corresponds to the 
question about annealed topological disorder, as it relates 
a n (t) to the probability of getting the topological state 
n. By contrast, the latter view on the same quantities 
a n (t) corresponds closer to the idea of quenched topo- 
logical disorder, as it reflects on the physical property of 
the polymer with given n. Another such physical quan- 
tity is the distance between the obstacle and an arbitrary 
point on the polymer. The exact expression for the root- 
mean-square of such distance is given by formula 14811 . 
The interesting aspect of this result is that this distance 
remains of the order of a\fi and only quite modestly de- 
pends on n, changing from approximately 0A96a\/i when 
n = 1 to 0.408a\/i when n — > oo. The fact that this dis- 
tance decreases with growing "topological complexity" n 
is not surprising, but the fact that it changes only slightly 
is interesting. One could have thought that the polymer 
would consist of n roughly similar loops, leading to the 
typical size of ay/t/n. Our result, therefore, suggests 
that even at very large n there remains one big loop, 
with the length of order t, while all other loops are tight 
and small. This is reminiscent of knot tightening recently 
discussed by Kardar and his co-workers |35| . 

In section lY Dl we make a brief comment on the elastic 
forces developing in the polymer ring either pushed too 
close to the obstacle or pulled too far away from it. 

Finally, in the section IV El we consider polymer ring 
entangled with a finite size obstacle, and show that in this 
case the distribution over the linking number n decays 
exponentially at large n, and the characteristic n is about 
ln(ta 2 /6 2 ). 



III. GREEN FUNCTION FORMULATION 

A. Point-like obstacle 

Consider a Gaussian polymer in 2D or, equivalently, 
a random walk in 2D. Suppose first that the obstacle is 
point-like, positioned at O, the coordinate center. The 
statistics of trajectories is fully described by the Green 



function, G 



which is the partition function (or 



statistical weight) of the chain having the monomer at 
r 1 and monomer t at r. The Green function satisfies the 
diffusion equation 
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where the notations are standard: A is the Laplace op- 
erator acting on r, a is the monomer size, t is polymer 
length ("time"). The notation a 2 /4 adopted here for the 
diffusion coefficient, which is in fact a 2 /2d, d being space 
dimension, is convenient because root-mean-squared end- 
to-end distance of the trajectory with no obstacles equals 
exactly a\ft. The Green function can be written in terms 
of the bi-linear expansion over the corresponding eigen- 
functions. Because our goal is to address the obstacle 
at O, we choose eigenfunctions with cylindrical symme- 
try. The ones with no singularity at O read J^nrje ^ , 
where J^(x) is Bessel function of the first kind, r and 9 
are polar coordinates corresponding to r, and —k 2 is the 
corresponding eigenvalue. Accordingly, we write 
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It is worth noting explicitly that only positive /i > 
contribute to this expansion, because J^{x) with negative 
index /i is singular at small x. 

In most cases in mathematical physics, the angular 9- 
dependence is 27r-periodic, meaning that 9 and 9 ± 2nn 
label one and the same place on the plane. This is not the 
case for the problem at hand. Indeed, G*(r, 0|r', 0) is the 
statistical weight of trajectories (polymer conformations) 
that start at a point some distance r away from the origin 
O and arrive after "time" t at another point some r' from 
O, where it is assumed that by the time t the trajectory 
has accumulated winding angle around O. Accordingly, 
for instance, 9 = means no turns around O, while 9 = 
2tt means one turn counterclockwise, 9 = — 2~k is one 
turn clockwise, etc. In other words, we should treat our 
plane as a Riemann surface, in which case 9 and 9 ± 2im 
correspond to different layers. 

Most immediately, this means that not only integer, 
but all positive values of /i must be included in the bilin- 
ear expansion Q. 

It turns out that integration over k can be explicitly 
performed; the derivation of the relevant so-called Weber 
integral [33] is provided in the Appendix^] The result 
reads: 
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where 7 M (x) is the modified Bessel function. 

It is instructive to re- write the latter formula by intro- 
ducing R - the distance between r and r* ': R = r — f f , or 
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The first factor in the Eq. |JBJ (in square brackets) is 
simply the Green function of an unrestricted polymer, 
or unrestricted random walk; in other words, it is the 
statistical weight of all conformations going from rior*. 
Therefore, W measures the fraction of trajectories with 
winding angle 9 on the way. 

Equations (|5I7|I were derived by Edwards [l4|. 



B. Finite size obstacle 

Consider now an obstacle having the shape of a disc 
with some finite radius b. Since the trajectory cannot 
make infinitely many turns around such obstacle, we ex- 
pect that the probability distribution for the number of 
turns should be completely different for this case as com- 
pared to the point-like obstacle. 

We use the same method as before. Eq. © still ap- 
plies, but as regards bi-liner expansion, Eq. Q , we have 
now different set of eigenfunctions - the ones which sat- 
isfy the boundary condition of being equal to zero at 
r = b. This boundary condition removes all trajecto- 
ries which cross the boundary, or, in other words, which 
enter the r < b region. The eigenfunctions, correspond- 
ing to the eigenvalue — k 2 , can be written in the form 
Z^inr, K&)e ve , where (see Appendix lB|l 
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Here Y^ (x) is Bessel function of the second kind (another 
frequently used notation for Y^{x) is Nn(x); we adopt 
here the notation used in Mathematica |3Jj|). A few notes 
about functions Z are provided in the Appendix [B| in- 
cluding the proof that the square root in the denominator 
makes them correctly normalized. Using Z, we write the 
Green function as a bi-linear expansion, like Eq. ifljl: 
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x Z^{Kr,Kb)Z^{K,r' , Kb) ndndfi . (9) 

Unfortunately, no known analog exists of the Weber in- 
tegral for the Z-functions, and so, unlike the 6 = case 
above, we were unable to find any way to simplify this 
by performing either of the two integrations. 

Addressing the sameproblem of winding around the 
disk, Rudnick and Hu |5| have already found the expres- 
sion for the Green function. Formula looks surpris- 
ingly different from the known result ■ In the Appendix 
|n]we show explicitly that these two results are equivalent. 



C. Winding inside the cavity 

Yet another interesting model is shown in Fig. [He). 
It is a random walk or linear polymer confined in a re- 
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stricted volume, say, inside the disc of some radius B. 
Then, absorbing boundary conditions should be imposed 
on this boundary. Assuming for simplicity that the ob- 
stacle is located in the center of the confinement disc, we 
obtain that Eq. 0} holds, except integration over k at 
every /j, must be replaced by the sum over the discrete 
spectrum of K n (fi) such that J^(/c„(/x).B) = 0. As usually, 
as t — ► oo we can resort to the ground state dominance 
principle, which means we can truncate the summation 
to one leading term: 

x J^r/B)J^r'/B)dn, (10) 
where £ M is the smallest root of the Bessel function J M (£). 

IV. WINDING ANGLE DISTRIBUTION: 
SPITZER LAW AND RELATED RESULTS 

A. Winding around a point (b — 0) 



The authors of the works UJ, la la IS tZl la la 113 , ex- 
amined the problem of winding angle distribution in the 
following formulation. Suppose the walker starts some 
distance r from the origin, and suppose we are interested 
in the winding angle distribution irrespective of r' , the 
distance from the origin to the trajectory end. Formally, 
such probability distribution is obtained via suitable in- 
tegration of the Green function over r': 
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In the Appendix [D] we show how to use the Weber inte- 
gral to follow this path. 

Unfortunately, in some other cases considered below, 
such as winding around a non-zero size disc (6^0), we 
don't have the advantage of the Weber integral simplifi- 
cation from Eq. to Eq. JSJ, which makes the explicit 
integration of the Green function over r' difficult. Be- 
sides, for polymer applications it is natural to keep track 
of the end position as long as possible. This is why it 
is useful to see how we can re-derive the Spitzer law Q 
directly from Eq. 0}, not resorting to Eq. (JSJ. This is 
what we shall do now. 

We note that integration over k in Eq. Q is effectively 
truncated at k 2 < A/ta 2 . When t is large enough, this 
leads to both kt and nr' being small. Then the Bessel 
function can be replaced by the first term of its expan- 
sion, Jfj,(£) — r(i+^) (l) ■ After that, the integration 
over k is easily performed, yielding 
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Assuming rr' /ta 2 <C 1 (see the discussion a few lines 
below), we see that the integral over fi is dominated by 



small fi in which area we can set T (1 + fx) ~ 1. In this 
approximation, the integration over \i is elementary, and 
results in 
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This is the Cauchy distribution for the winding angle 
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This result is similar, but not identical to the Spitzer 
formula Q). The difference is in the definition of the 
scaling variable x: formula 114(1. unlike (IJ, keeps track 
of the coordinates r and r' of the trajectory ends. As we 
have already mentioned, this will be useful for polymer 
applications. Note, however, that we cannot integrate 
over r' as in Eq. i|ll|) . because formula l|14|) was derived 
under the assumption that r' is not too large. 

How then can we recover the Spitzer law Q from Eq. 
(I14JI ? What we should do is to note that one trajectory 
end is fixed at the distance independent of t, while the 
other is free, meaning that r ~ a and r' ~ a\ft. Then, 
we have for the scaling quantity x in formula (|14|) x = 
9 /in (ta 2 /rr') ~ 20/ kit, wkich is indeed exactly the 
same as in Eq QJ. 

Other interesting extremes are as follows: 

• If r ~ a and r ~ at, then the "width" of the dis- 
tribution gets very small. This is the closest ap- 
proximation Gaussian model can provide for the 
idea that fully stretched polymer does not have any 
freedom to wind around the obstacle. Of course, a 
Gaussian polymer cannot be fully stretched, this 
is why, say, {9 2 }, remains divergent even when the 
"width" goes to zero. 

• A similar situation is realized when r ~ r' ~ a\ft: 
winding is suppressed when the obstacle is removed 
to the periphery of the coil. Note that Eq. 1(14(1 
should not be used at larger r, when kt is not small 
and the Bessel function cannot be expanded. 

• If both r ~ a and r' ~ a, then the result is only 
different from Eq. P by a factor of 2 in the defi- 
nition of x; in this = 9 /hit. That means, 
fixing both ends and not allowing them to wander 
freely reduces the "width" by half. 



B. Winding around a disc (6 > 0) 

For winding around a disc of finite radius 6, we can use 
the same method. When t is large enough, integration 
over k in Eq. @ is dominated by small k. Accordingly, 
we can resort to the small k expansion of Z M (Kr, Kb) (see 
Eq. ((B19() and the discussion in the Appendix IB} : 
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Accordingly, the K-dependent factor in the Green function JUJ) can be presented in the form e g ( K \ where 
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Provided that /i < 1 (which is justified a few lines below), it is not difficult to establish that g(n) has a minimum, 
which dominates integration over n at large t. Straightforward differentiation yields for the the corresponding k the 
condition 



«V = 4,- (¥)- 2 "r»( 1 + fl ,-( ¥ )-r»( 1 - fl , 



(f ) 2M T 2 (1 + /i) + (f )^T 2 (1 - ,i) - 27r/XCOt7T/Z 



(17) 



This equation has just one solution which at large i cor- 
responds to small k. More accurately, the solution reads 
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1 when /iln(ia 2 /6 2 ) < 1 
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As it turns out, the integral over fj, is dominated by /j <C 
1 /in (ta 2 /b 2 ) , so only the lower line of the Eq. (|18fl is 
relevant. For small /i, the expression for Eq. (|15f) can 
be further simplified: 
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Then, replacing e 9 ^ K ' with its maximal value, we arrive 
at the following expression for the Green function: 
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neglected. That means, the convergence of the inte- 
gral H2()(l is controlled by the sinh in the denominator, 
which effectively truncates integration at (x smaller than 
1/ln (ta 2 /b 2 ). This is very good news. First of all, since 
l/ln(ta 2 /6 2 ) <C 1, this justifies the small /i simplifica- 
tion performed in formula i|19|) . Second of all, this also 
justifies the use of the lower line in the expression (|18ll for 
the saddle point. Third, since only small /i contribute, 
the validity condition for the expansion of Bessel func- 
tions in the first step of Eq. H15fl , which generally reads 
{kt) 2 < l + /i, can be simplified to r 2 <C a 2 t (and simi- 
larly r' 2 < a 2 t). 

Thus, all approximations leading to the expression H20(l 
are self-consistent. The only task left is to evaluate the 
integral (|2U[) . This task gets easier if we use the notations 
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Then, formula l|20|l is transformed into the following ex- 
pression for the probability distribution of the winding 
angle 9, or, better, of the scaling variable x, at fixed r 



plus some logarithmic corrections. In A we accumulated 
all the uninteresting constant prefactors, which do not 
depend on 9, r, r', and b. 

Now, considering this integral over fi, we have to jus- 
tify all the assumptions and approximations which we 
made on the way. First and foremost, the assumption 
that fj, is small is justified by the rapid convergence of 
the integral (|20|l . Indeed, at large fi all three sinh's 
can be replaced by positive exponentials, leaving us with 



exp 



'In- 



-lnl-ln^ 



Since 



< ta 2 , the 



latter two logarithms in the round brackets should be 



W(9) 



2aa' 



+OC 



sinh (£a/ir) sinh (£a /tt) 

'■ — 2 1 

sinh £ 



(22) 

where we have re-introduced the normalization factor, 
such that W(9) dx = 1. This integral can be re- 
duced to the infinite sum of residues corresponding to the 
poles along the imaginary axes on the complex £-plane. 
In turn, the resulting sum (which is the combination of 
several geometric series) is easy for Mathematica [33| , but 
can be also computed by hand. One way or the other, 
here is the result: 
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W{6) 



ttx sinh tt x sin a sin a' + a sin a' [cos a' — cos a cosh ttx] + a' sin a [cos a — cos a' cosh ', 



2aa' 



[cosh7ra; — cos (a — a')] x [cosh7rx — cos (a + a')] 



(23) 



The result (|23jl is unfortunately quite cumbersome, al- 
though it is symmetric and in some ways quite nice. Its 
beauty is revealed by consideration of various limits. As 
we learned in the case of point-like obstacle, the most 
interesting limit is when chain end is free, meaning that 
r' ~ ayi. Then, provided only that ta 2 ^> b 2 - which is 
necessary, as the walker must have traveled much farther 
than the obstacle size 6, we get a' ~ tt. In this case, we 
get 



W{9) = - 



sin a. I a 



vr/4 



2 cosh ttx + cos a 



cosh (irx/2) 



(24) 



where in the latter transformation we also noted that as 
the trajectory starting point is fixed, r is independent 
of t, or a — > at large t. Thus, we recover formula 
(J2J) • Importantly, the definition of scaling variable x (|21|) 
becomes identical to that in (|2J), again under the same 
condition ta 2 » b 2 . 

As in the case of point-like obstacle, other interesting 
extremes are as follows: 

• Both r ~ ay/t and r' ~ ayt. In fact, this case is on 
the border of applicability of our approximations, 
but qualitatively the result holds. Indeed, W(9) 
becomes very narrow, and approaches 8(x). This 
means, no turns are possible around the obstacle 
which is away from the random walk trajectory. 

• Another case, and also a border case in terms 
of applicability of our approximations, is r ~ at 
or r' ~ at, implying an exponentially improba- 
ble straight trajectory. The distribution is again 
sharply localized at small x. 

• Both r and r' are independent of t, meaning that 
both a ~ and a' ~ 0. This case is safely within 
the limits of applicability. Then, 



2(1 — cosh-; 



„,,.. TT TTX Sillh TTX 

W{ ~ 9) = ? n v, — W 

*> ( 1 — cosh ttx) 



(25) 



As in the b = case, this distribution, as one can 
easily check, is exactly two times more narrow than 
that of (J2J. In this case, unlike b — 0, this statement 
can be formalized by looking at the second moment 
of the distributions (J5J) and l|25() . which (in terms 
of x) turns out equal 1/3 and 2/3, respectively. 



6^0 limit: applicability conditions of the 
Spitzer formula 



variance at 6 = 0. These equations leave it unclear what 
happens when the obstacle gets smaller and smaller, or 
when b decreases and approaches 0. It is instructive and 
interesting to use formula 123(1 to sec what really happens 
when b — > 0. 

The important part of our analysis here is to realize 
that so far we have been using several differently defined 
scaling variables x: see Eqs. (Q), ©, JT3J|, l|2*Tjl. So far, 
it was (hopefully) clear from the context in every place 
which x we have in mind. Now, when we examine the 
b — * limit, we shall face several of these different x 
simultaneously, so we must be certain as to which x is 
which. For the rest of this section, we adopt the notation 
in which each x is labeled with the number of the defining 
equation: a f^t (note, that aj^jis exactly the same 
as aj^). In particular, x m formula l|23(l is, of 

course, aj^fl 

Speaking of different definitions of x, we should realize 
that so far we have been presenting probability distribu- 
tions W{9) normalized with respect to integration over 
the corresponding x. For our purposes now, it is more 
convenient to use the normalization condition with re- 
spect to angle 9: f™^ W(9)d9 = 1. For the formula P5fl. 
this means the factor 2/ dn(ta 2 /b 2 ) should be incorpo- 
rated; we do not re-write the formula for brevity. 

The main reason why the difference between various 
x was unimportant so far is that at t — > oo all defini- 
tions converge to the same: aff^l^ 21^— > am= 29/ Int. 
However, when b — > 0, there appears a very broad inter- 
mediate range of times t such that although ta 2 ^> b 2 (the 
trajectory is long enough to wind around the obstacle), 
but tb 2 < a 2 : 



a 2 jb 2 » * » b 2 /a 2 



(26) 



This is the range which we must examine. In this range, 
to the leading approximation, a f^l does not depend on 
time: 



a EJ= 



29 



ln(o/6) 



(27) 



Furthermore, there is a broad range of winding angle 9 
in which a|2T] i s small. 

Now, we should look at the quantities a and a' . When 
6 — > 0, both of them turn out to be slightly below ir. For 
instance, 



2tt ln ; 



lnf 



In* 



b 



ilnt 



(28) 



where 



According to the Eqs. @J and J2J, winding angle distri- 
bution has finite variance at b ^ and diverging infinite 



2 lnf 



< 1 



(29) 
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Similarly, a' = tt — 5', with similarly defined 5' -C 1. 

Thus, we can simplify formula (|23[1 resorting to ex- 
pansion of both numerator and denominator over the 
powers of a|2"Tl <^> and • I n fact; as we see from 
Eqs. H27I28|) all these expansions are ones over the in- 



verse powers of In (a/ b), and we keep the leading terms 
only. Incorporating, as explained above, the factor 
2/ ln(ta 2 /b 2 ) ~ l/ln(a/6) to establish the normalization 
f-oo w ( d ) de = 1, we finally get 



W{6) 



^(5- S'f (S + 5>f + ^[(5 + Sf + (5- Sf\ 



hl !l) l n *4#2 



In- 



In T - 



V rr / 



9 2 +9 4 



(inK) 

\ rr' J 



9 2 



(30) 



For the two latter steps, we have plugged in the explicit 
expressions for S, 5' lj2l31) . and ljT7|) . and then ne- 
glected the d — 6' ~ ln(r'/r) terms (which we have un- 
derbraced in the intermediate formula). The result is 
exactly the same as formula Ijl4l) (except it is normalized 
with respect to integration over 9). 

From our analysis, we can now understand the cross- 
over between Eqs. J5J and Specifically, the Spitzer 
formula and its generalization (|14(l apply as long as 
two conditions are met: t <C a 2 /b 2 and 9 <C ln(a/b). At 
longer times and/or larger angles, the exponential tail of 
the distribution takes over. For instance, when wc say the 
(9 2 ) diverges for winding around a very small obstacle, 
this really means (9 2 ) ~ (ln(a/6)) 2 . 



D. Winding inside a cavity 

We start from Eq. (|10|) . It is not difficult to realize 
that £ M (the smallest zero of J^i^)) increases with /z. 
Therefore, when t gets large, the integration over /i is 
dominated by small //, as in all previous cases. At small 
fi, £p is a smooth non-singular function, we can linearize 
it: ^ ~ £o+M£o- Numerically, Co ~ 2.405 and £' w 1.543. 
To the same approximation, J^^r/B) ~ Jo(Co r / B). 
Therefore, evaluation of the integral in Eq. I|10|) becomes 
trivial, and the result reads 



G 



r',0 




t 



Hf) J <^) wm ' (31) 

where the probability distribution of the winding angle 
is given by 

ta 2 ^' /2B 2 



1 



W{6) = - 



(32) 



n(ta 2 ^' /2B 2 r + 9 2 

The decoupling of the ends r and r' in formula (|31|l is 
not surprising, this is the property of random walk locked 



in a restricted volume, and it is due to the fact that 
correlations are broken every time that the trajectory is 
reflected from the cavity border. As regards probability 
distribution of winding angle, it is once again the Cauchy 
distribution, however, the scaling variable involves 9 j \ft 
instead of (9/lnt for the random walk in an unrestricted 
space. This is also because correlations are broken every 
time that the trajectory hits the border. One can say 
that pieces of random walk with length about (B/a) 2 
act independently of each other. 

This gives rise to the following simple scaling argument 
providing an insight into the result (|31(l . The winding an- 
gle distribution for every "blob" of the length ~ (B/a) 2 
is given by the Spitzer formula Q , with the replacement 
t — * (B/a) 2 . Now, we have tj (B/a) 2 of such blobs. Since 
blobs are independent, the probability distribution of the 
sum of all winding angles of all blobs is given as a convo- 
lution. In other words, Fourier transform of the Spitzer 
distribution for one blob, which is e - ' M ' ln ( B / a ) , must be 
taken to the power t/(B/a) 2 . Apart from logarithmic 
corrections, this returns the result 131|) . 

Thus, the reason why (9 2 ) diverges for the polymer 
inside the cavity is because every blob can make many 
turns around the point-like obstacle on a small scale, be- 
fore ever hitting the border of the cavity. 



V. RING POLYMER: 
EDWARDS-PRAGER-FRISCH MODEL 

We now want to make one step closer to the attempt 
of gaining insight into the properties of closed ring poly- 
mers. One way in this direction would be to say that a 
ring polymer is the random walk trajectory whose end 
points happen to coincide, namely r = r and 9 = 2nn, 
where integer n (positive, negative, or zero) is the linking 
number (number of turns) . Our results 114I23I25I32|I are 
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suitable for this, and we shall do it. It turns out also 
useful, however, to derive some additional results inde- 
pendently. In particular, some of the results below are 
exact (not asymptotically exact, but just exact). 



Point-like obstacle: a ring with one monomer 
anchored 



Thus, we return to Eqs. (|6I7[1 . and use them this time 
to write down the statistical weight of the ring poly- 
mer conformations with linking number n and with one 
monomer fixed at the distance r from O: 



G.A.r.h = < ( ^ 



r, 2irn 
t 



-4^^ i ( 33 ) 
7ra z r 



f°° 2r 2 
2e~« / cos(2?m/x) 1^) d/j, ; £ = — . 
Jo aH 



In this formula, G n (r, t) is the statistical weight of the 
ring with n turns, while the prefactor l/ira 2 t is the sta- 
tistical weight of a ring with no topological constraints 
|4l| . Therefore, W n {£) is the probability that polymer 
ring fixed at one point r makes n turns around the ob- 
stacle. In the Appendix [E] we check explicitly that W n 
satisfies the normalization condition as a probability. 

Similarly to what we did before, we can address the 
case ta 2 r 2 , or £ -C 1. In this case, we truncate the 
small £ expansion of 1^(0 — (£/2) /1 /r(l + /_t), replace 
r(l + fj.) ~ 1 (compare Eq. I|D1||1. and then obtain 



W n (0^ 2(1-0 



ln(2/0 



(ln(2/0) +47r 2 n 2 



^0. (34) 



Of course, this is nothing else but the "discrete" version 
of the Spitzer distribution. However, merely taking 9 = 
2irn in Eq. Q14f) is not enough, as the normalization factor 
in ltTl|) corresponds to /_ °° W(6)d6 = 1, while in Eq. 

(El it corresponds to W n(Q = L 

For a polymer, it makes perfect sense to examine also 
the opposite extreme, £ 3> 1. The corresponding asymp- 
totics are easy to derive from the somewhat simplified 
expression for W„(£). 

We can afford further simplification of the expression 
Eq. ljo"4"|l for W n (£) resorting to the following integral 
representation of the modified Bessel function 



U0 = T- r ^ cos/3 cos(^)d/3- 

27r J -ix 



sin(/i7r) 



(35) 



which is generally valid at |Arg£| < ir/2 and 5R/i > 0. 
Both of these conditions are met in our case. Substitut- 
ing this into the Eq. (|34|) . one can easily perform the 
integration over fj, yielding 



W n {0 = A„ 



du e -«( 1+coshu ) x 



(36) 



2n- 1 



2n + 1 



v 2 + 7T 2 (2n - i f u 2 +tt 2 (2n+ l) 2 



where A„ is the Kronecker symbol (1 for n — and 
otherwise) . The latter result for n = is worth re- writing 
separately: 

r°° p -?(l+coshu) 

W (0 = l-2 -^^du. (37) 



l 
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FIG. 2: W n (£,) is the probability to form a link of order 
n with the point-like obstacle provided that one of the chain 
points is fixed at the distance r away from the obstacle, where 
£ = 2r 2 /a 2 t, t being the chain length. The plots present the 
result of numerical integration based on formula 1361 . The 
plot of Wo(£) is presented in the main figure. W„(£) with 
7i > would not be seen well in this scale. For both Wi(£) 
and W2(£), we show the inset, each presenting the vicinity of 
the maximum; the corresponding places on the main figure are 
shown by tiny dark rectangles. However small may seem every 
particular W n (Q, it should be born in mind that together they 
sum up to 1 — Ho(£). 



Equations (|36I37|I are convenient enough to address 
the £ 3> 1 extreme. Indeed, when £ is large, the integral 
converges at small u, which allows us to neglect u ev- 
erywhere except in the exponential factor, where we can 
also truncate the coshu ~ 1 + u 2 /2. This yields: 



4n 2 -l 



v 7 ? 

1- V^Ar 



3e_ 



n = 



(38) 



Thus, comparing and (jSEl all W„(£) start at at £ = 
and grow very rapidly at small £. At n ~ 0, Ho(0 keeps 
increasing monotonically with £, and Wo(£) approaches 
the saturation level of 1, while all W n (£,) with n ^ 
decrease and rapidly die away at large £. Obviously, each 
of them goes through a maximum. It is not difficult 
to establish that the maximum of W n {^) corresponds to 
£ ~ 1/ cosh(7T\/4?T. 2 — 1) which at large n corresponds to 



-2-Kn 



This is consistent with the fact that small £ 



asymptotics Eq. (P3} is valid at £ < e~ 2?m . 

Eq. I|36|) allows also straightforward numerical integra- 
tion which results in the plots shown in the Fig. [21 
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B. Ring polymer entangled with a point 

For the ring polymer, it is not very natural to consider 
one monomer being fixed at r; all monomers of a ring 
are equivalent. Accordingly, it is natural to define the 
quantity 



/>oo 7ra 2 t f°° 

Vn(t) = / W n ^)27rrdr = — / W r , 
Jo 1 Jo 



(39) 



What is <7„(i)? This quantity has the units of surface area 
and can be interpreted in the following way. Suppose 
a ring polymer moves freely on the plane within some 
large area A (much larger than the polymer size, so the 
polymer is not restricted in terms of its conformation). 
Consider one particular conformation of our polymer and 
then choose a random point O within A. Then cr n /A is 
the probability that polymer makes an n-fold link around 
O. We expect physically that er should be large, almost 
as large as A. This is also seen directly from the Eq. 
(|37|l : when integrated over the whole area A, the first 
term (unity) yields just A. This is because when O is 
outside the coil, there may not be any topological links. 
If and only if the random point O is located within the 
polymer coil can there be any topological link. Therefore, 
(To should be less than A by a quantity of the order of 
the coil gyration radius squared, which is of the order ta 2 . 
On the other hand, o~ n with n ^ should be themselves 
of order ta 2 , or even smaller. 

Another way of understanding of a n (t) is this. Con- 
sider one particular conformation of a ring and consider 
the set of points S„ such that the polymer makes linking 
with n turns around every point of S n . Then, cr n (t) is 
the surface area, or the measure, associated with the set 

Trying to compute o~ n (t), we can resort to either of 
the expressions (|54T) or Let us first explore the first 
possibility: 



roo poo 
,2 / -£ 



o- n (t) = nta 2 e -6 / cos^tto/x)!^)^*; • (40) 
Jo Jo 

Here, we face a difficulty, because the integral 

/>oo 

/ e"%(£K (41) 
Jo 

diverges at large £. What is the physical meaning of this? 
Of course, this is because <Jo(t) is close to A, or, in other 
words, it is divergent unless we take into account over- 
all volume restriction. We conclude, therefore, that the 
integral l|41|l diverges for a good reason: this is because 
unlinked polymer is free to move away from the obstacle, 
making ctq a s large as (almost) A. 

This hints on the way to circumvent the problem. Let 
us assume that the polymer is attached to the point O 
by a very weak spring. Since such polymer does not 
move away even when there are no topological links, we 



expect that even o~q should remain finite, independent of 
A. Indeed, instead of (|41|l we have now 



fV Q %(£K = 
Jo 



1 



\Ja 2 - 1 (a + \/a 2 - l) 



(42) 



which converges at any a > 1 ; here a — 1 is the effective 
spring constant. Of course, we will take a — > 1 at the 
end. Performing the remaining integration over fi, we 
arrive at 



a n (t, a) 



nta 2 In (a + \J a 2 — l) 



y/a^T (27rn) 2 + (\n (a + Va 2 ^!))' 



(43) 

As expected, the a — ► 1 limit can now be performed with 
no difficulties at every n^O, yielding finally 



. . ta 2 1 

0- n (t) = j i 11 T 

47T 



Accordingly, 



a (t) =A-J2<rn(t)=A- ^Ua 2 

n/0 



(44) 



(45) 



That result exactly can be also obtained plugging Eq. 
Ijflfijl into the Eq. . although, somewhat surprisingly, 
calculations are more involved along this route. 

We would like to remind to our reader once again, that 
o~ n (t)/A is the probability to have linking number n for 
the polymer of the length t. As the probability distribu- 
tion, o~ n has the peculiarity that all its moments obviously 
diverge, even just the average linking number is infinite. 
It is not difficult to trace this back to the fact that in- 
finitely flexible polymer, as represented by the Brownian 
random walk trajectory, can make infinitely many turns 
around a point-like obstacle. We shall address this fur- 
ther later. 

It is worth emphasizing that the results (|44I45|I are 
exact, their validity does not require even that t is large 
- they are exact at any t. 



C. How far is the ring from the point-like obstacle? 

One more interesting quantity to look at is (r 2 ): the 
mean squared distance of one particular point on the ring 
to the obstacle, O. To determine the probability distribu- 
tion for r, we note that <r n (t) plays the role of a partition 
function. The probability density for r reads 



G 


( 0,r 

{ o 


27m, r 
t 


) 




' 0,r 



2nn, r \ 
t ) 


2-irrdr 



0"n(i) 



(46) 



When computing (r 2 ) from this, formula (|43|l comes in 
handy, as (r 2 ) is basically the derivative of a n (t, a) with 
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respect to a at a = 1: 



E. Ring polymer entangled with a finite size disc 



(r 2 ) = 



7r ( a f\ 

it) / w n (0^ 



a 2 t d<j n (t,a) 



2<r„ (t) da 
Straightforward calculation yields 



aH 



< r > - - 



1 



2?m 2 



(47) 



(48) 



The result is interesting. Surprisingly, it goes to a finite 
constant proportional to the unperturbed coil size ta 2 in 
the limit of very strong linking, n — > oo. This should 
be understood by noting that even very large number of 
turns will be produced by a short piece of a trajectory, 
leaving a long part, of the order t, unentangled, with the 
size of order a 2 t. 

This is reminiscent of the recent findings by Kardar 
and his co-workers |35| in which they claim that in many 
cases real knots in three dimensional polymeric loops are 
entropically dominated by conformations with the knot 
tightened in a short piece of polymer, and with the rest 
of the polymer fluctuating freely with no knots. 



D. Force 

When windin g m odel was first introduced in the poly- 
mer physics |l4lll5llT^| . it was done mostly in connection 
with problems related to the rubber elasticity. Accord- 
ingly, elastic force, or force-extension curve, was the pri- 
mary subject of interest. In case of DNA, such a curve 
can be also measured using some sort of a single molecule 
technique |3fil |. Although both in rubbers and in DNAs 
real forces have both entropic and enthalpic contribu- 
tions, in the winding model the force has purely entropic 
nature and, therefore, it is proportional to ksT in stan- 
dard notations, where ks is Boltzmann constant and T 
is absolute temperature. In our notations, the force /„ 
which should be applied to the polymer to keep one of 
its links a certain distance r from the obstacle O at the 
fixed topological invariant n is given by 

fn dlnW n 4r dW n 

k B T dr ta 2 W n d£ ' [ ' 

We have not found any simple closed expression for the 
force, in this sense we make really no progress on this 
point compared to the papers [lj, [jja, [l6| • Nevertheless, 
qualitatively, one glance at the Figure |21 is sufficient to 
realize that the force /o of an unentangled ring is always 
positive. This is obviously because this ring is topolog- 
ically repelling the obstacle. On the other hand, when 
n ^ 0, the force is positive, corresponds to repulsion, only 
when r (or £) is small enough. At larger r (or £), the force 
flips sign and becomes negative, which obviously corre- 
sponds to the elastic stress caused in the polymer ring by 
an attempt to pull it away from the obstacle with which 
the ring is entangled. 



For the obstacle of finite radius 6, we were unable to 
obtain exact answers similar to Eqs. I|44|) or 1)48(1. All we 
can do for this case is to resort to the asymptotic calcu- 
lations. One of the advantages of the finite size obstacle 
model is that it allows to examine both asymptotics, we 
call them loose entanglement and tight entanglement, re- 
spectively. The former regime is realized when the size 
of the obstacle b is smaller than typical polymer coil di- 
mension at 1 / 2 , and, moreover, when minimal length nec- 
essary to make n turns, 2nbn, is still small compared to 
at 1 / 2 : nb <C at 1 / 2 . In this case, calculations are similar 
to those of Section TlV Bl In the opposite extreme, when 
nb at 1 / 2 , polymer has to be significantly stretched out 
to make all n turns. This corresponds to the far tail of 
winding angle distribution, which is usually not exam- 
ined and which we did not consider in Section llV Bl 



1. Loose entanglement 

All we can do for this case is to resort to the asymptotic 
calculations similar to those of Section liVBI In fact, the 
calculations are almost identical, and at the end they 
return essentially the result with the only difference 
in the normalization factor. Specifically, the probability 
to have linking number n is proportional to 

jjt nx n sinh7nz;„ + 2 (1 - cosh7ra;„) 

W n cx -g . (50) 

( 1 — cosn irx n ) 

where the omitted normalization factor must be defined 
such that XmL-oo = 1) an d where 



Ann 



In 



ta 2 
b 2 



(51) 



Similarly, although we cannot find the exact expression 
for the value of o~ n (t), but the estimate reads o~ n (t) ~ 

Thus, quantities such as W n and a n (t) decay expo- 
nentially at very large n, and the characteristic n where 
exponential decay starts is about ln(ta 2 /6 2 ). This latter 
quantity estimates also the characteristic linking num- 
ber in another sense, defined as the root mean squared, 
\J (n 2 ). This is an interesting and somewhat unexpected 
result. Indeed, one could have naively expected that the 
characteristic value of n should be proportional to ta 2 jb 2 . 
Indeed, we expect that one turn around the obstacle 
should be similar to walking a distance about 2-Kb ~ 6, 
which requires a time r ~ b 2 /a 2 , implying the number 
of turns to be about t/r. Instead, we are getting some- 
thing like ln(t/r). This happens because a large portion 
of the chain length deviates much further away from the 
obstacle than &, and it makes turns around a much larger 
circumference. This once again suggests that knot tight- 
ening [35| occurs even in this case of a disc-like obstacle 
with excluded volume. 
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Similar to our discussion in section IIV CI we can un- 
derstand what happens when b — > 0. In this case, 
there appears a wide interval of polymer chain lengths 
b 2 /a 2 <C t <C a 2 /b 2 in which, say, <J n {t) decays only as 
1/n 2 in the wide interval of n, up to a large value of n of 
about ln(a/6). 



2. Tight entanglement 




The tight entanglement regime is realized when n 2 b 2 
ta 2 . In this case, polymer barely has enough length to 
make n turns around the obstacle. Obviously, the dom- 
inant polymer conformations are those tightly wound 
around the obstacle. This regime is similar to ray optics 
|38j . Indeed, if one searches for the solution of diffusion 
equation in the form G = exp [s(r, t)] and assumes 
that s is (in a proper sense) a slowly changing function, 
then the so-called eikonal equation for s results: 



ds 
dt 



G 



Vs 



(52) 



For the system at hand, namely a polymer with winding 
angle 8 around the obstacle of radius b, this equation 
allows for the exact solution: 



2a 2 t 



(53) 



where L is the shortest distance between fixed ends con- 
sistent with the topological constraint (that is, with the 
given winding angle): 



L 



$-§') + \/r 2 - b 2 + vV 



b 2 



(54) 



Here, i? and are determined by the conditions cos# = 
b/r and cost?' = b/r'. Both these conditions and the 
solution itself are quite easy to establish based on the 
geometry presented in the figure^ It is also easy to check 
by direct differentiation that formulas (|53I54|I present an 
exact solution of the eikonal equation Q52[l. 



F. Ring polymer inside a cavity 

Our discussion in the previous section is additionally 
illuminated by the problem of a ring polymer entangled 
with an obstacle while confined in a cavity of the radius 
B. In this case, the result directly applies, apart 
from the replacement 9 — > 2im, and the proper normal- 
ization factor. We see that in this case the characteristic 
value of n is proportional to ta 2 /B 2 . This number can 
be understood by saying that polymer is confined to a 
tube of the width D = B and length L — nB, where the 
typical n must be determined such that confinement en- 
tropy is similarly contributed by chain squeezing across 
the tube and stretching along the tube (compare simi- 
lar arguments in [39)). The confinement entropy is well 



FIG. 3: Tight entanglement, or "ray optics" limit. In this 
figure, for the ease of drawing, we assume that the polymer 
makes just a little more than one turn around the obstacle, 
while its ends are fixed at the given points. Distance to one 
end r and the corresponding angle ■& are shown in the figure; 
similar distance r' and angle are not shown to simplify the 
figure. 



known ^£j> it is proportional to ta 2 / D 2 + L 2 /ta 2 , where 
the two terms correspond to the two factors just men- 
tioned - chain squeezing across the tube and stretching 
along the tube, respectively. Equating the two terms, we 
arrive at n ~ ta 2 /B 2 , as expected. 



VI. CONCLUDING REMARKS 

We have focused on a s peci al application, a toy 
Edwards-Prager-Frisch model |l4Lll5l |. and its generaliza- 
tion. It can be viewed as a model of an unsymmetrical 
infinite catenane. This catenane is formed from a ran- 
dom walk of t steps and "entwined" with infinite, rigid, 
closed structure composed of two straight legs, which are 
separated at least by a distance larger than ta, meeting 
at infinity. As such it allows one to extrapolate a limiting 
probability of catenation by a closed random walk, which 
is consistent with earlier estimates. 

This model is widely recognized as the simplest play- 
ground for "statistical mechanics with topological con- 
straints" Unfortunately, no simple notable result 
had previously come out of this model studies - except 
the very fact that it is "exactly solvable." We were lucky 
to find a couple of such simple results. First, looking at 
the entanglement as an element of "annealed" disorder, 
we found that the area associated with all points around 
which closed random walk makes exactly n turns is equal 
to (R 2 ) /Ann 2 , where (R 2 ) is the mean-square end-to-end 
distance of the linear walk of the same length. Second, 
looking at the entanglement as an element of "quenched" 
disorder, we found that the mean-squared distance be- 
tween an obstacle an an arbitrary monomer of an n times 
entangled ring is equal to (1/6) (R 2 ) [l + 3/27rn 2 ] . Both 
results are exact. We have also found that the entan- 
glement of a very long polymer is very uneven, in the 
sense that it tends to segregate into one very long loop, 
almost as long as the entire polymer, and a number of 
much shorter loops. 

The generalization of Edwards-Prager-Frisch model, in 
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which an obstacle is not a point, but a disc of finite radius 
b, as far as we could tell, does not allow for an exact 
solution in any useful closed form. However, we were able 
to show that for annealed loop the mean-squared linking 
number turns out to be finite, of order (n 2 ) ~ ln(a/6). As 
b increases, the entanglement becomes increasingly tight, 
resulting in all loops being of comparable length. 

Our deliberations clearly reflect three aspects: The 
first is the fruitful nature of Spitzer's [1| original insight 
for both further theory and its applications. The second 
is the value of self-consistent approximations in the field 
when direct exact calculations are precluded by mathe- 
matical difficulty. And, last but not least, the third is the 
usefulness of complete analysis of those simplified models 
which do allow for mathematically exact solutions. 
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which is most frequently encountered as an integral rep- 
resentation of 1 
lation leads to 



resentation of the modified Bessel function J M . This re- 



J^y)I,,{2a 2 xy) x e a v ydy 



2a 2 



(A4) 



Changing variables, we finally obtain the two formula- 
tions of the Weber integral: 



x e a y ydy 



e 4 



2_ K 2 
2 



2a 2 

. 2„ 2 



-Ju 



" \2a 2 

J^y)J^'y) x e- a ' y 'ydy 



2a 2 "MS' - (A5) 



The latter formula is what needs to be used to go from 
Eq. m to Eq. ©. 



APPENDIX B: SOME PROPERTIES OF THE 
FUNCTIONS Z 



APPENDIX A: WEBER INTEGRAL 



1. Orthogonality 



First let us prove an auxiliary relation: 

-c?{f-f r^jdV = ^e-^ 4a2 Mf) , (Al) 



where ip K (r) is an eigenfunction of the Laplacian operator 
corresponding to the eigenvalue —k 2 : Aip — —K 2 ip(r). 

To see this, we note that (a 2 /ir)e~ a ( r ~ r ) is the Green 
function of the diffusion equation in which a 2 plays the 
role of 1/Dt, with D and t being diffusion coefficient 
and time, respectively. Therefore, this exponent can be 
written in terms of a bilinear expansion 



a/- -''i 2 



k 



ipk{r)?Pk(r )e 



(A2) 



where the summation runs over the entire spectrum of 
the Laplacian operator. By making a dot-product of both 
sides with ip K) we arrive at the result IjAljl . 

Let us now use the formula (|A1|) choosing ip K = 
J n (nr)e m6 . At any n, this is indeed one of the eigen- 
functions of the Laplacian operator, corresponding to the 
eigenvalue — k 2 . Now, we make use of the following rela- 
tion 



/■27T 



OL COS <fi-\-lfA<fi 



2nl^(a 2 ) 



(A3) 



First of all, we want to prove here that the functions 
Z^kt, Kb), as defined by the Eq. ©, are orthogonal and 
normalized: 



1 

Z„(Kr, Kb)Zn(Kr, K'b)rdr — —S(k — k) (Bl) 

k 

1 

Z^kt, Kb)Z^(Kr' , k6)kc?k = — S(r — r ) (B2) 



Note that when b — > 0, we have Z^(Kr,Kb) ~ J^(kt) 
(because in this case Y^Kb) is negative and large in ab- 
solute value), so in this limit both equations (|B1|I and 
(IB2I) come back to the well known relation 133 



J ^(Kr) J ^{k 'r)rdr = —S(k — k ) . (B3) 



In what follows, we derive the relations (|B1IB2|) . For 
simplicity of notations, it is easier to compute the nor- 
malization of the functions 

U^Kr, Kb) = -J^Y^Kb) + J M (/t6)y M (/sr) , (B4) 

from which the properties of Z will follow automatically. 
Derivation consists of two parts, one of them is trivial, 
and the other is only slightly less trivial. 
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a. Trivial part 
To begin with, there is a useful general formula 



U 2 (nr, Kb) r dr 



9 2 
r _]L_Tj2 

2k 2 " 



r 

2^? 



dr 



(B5) 

which is valid for any solution of Bessel equation, that 
is, for any linear combination of J^(kt) and Y^(nr), in- 
cluding the Up. The derivation of this formula can be 
found in many places, for instance, |37j . One way is to 
take (rdlf^/dr) 2 and differentiate it over r. Remember- 
ing that satisfies Bessel equation, it is easy to find 



0_ 

Or 



( dUA 2 2 


— - K 2 


{ r dr) = 2r 





or 



(B6) 



from which Eq. (|B5|I follows automatically. 

We cannot directly apply this formula for the case of 
an infinite interval, because it yields the divergence (due 
to the r 2 U 2 term: decays only as ly/r at large r). 
Indeed, this is not surprising, since the answer (|B1|) con- 
tains a ^-function. Thus, what we shall do is to consider 
first the finite width ring b < r < B, with boundary con- 
dition U f2 \ r=B = 0. In the end, we shall send B — > oo. 

Assuming U^nr, nb)\ B , we get from <|B5|) 

B 



U 2 (nr, nb)rdr 



r 

2^2 



dr 



(B7) 



J b 



The derivative dU^/dr can be simplified, because it is 
related to the Wronskian of J^nr) and Y^nr), which 
is equal to 2/nnr. Taking into account the boundary 
condition 



Y^ K b)J^(KB) = Y„(/sB)J /t (/s6) , 



(B8) 



we obtain 



U 2 (nr, nb)rdr 



J 2 M) 

J 2 ( K B) 

Y 2 {nb) 
Y*(kB) 



1 



(B9) 



If we have two different values, k ^ k', both satisfying 
boundary condition (|B8() . then 

rB 

U^nr, K^U^n'r, n'b)rdr = , (BIO) 

as can be established either by proper integration by 
parts using the Bessel equation, or by direct reference to 
the fact that these t/'s are the eigenfunctions of a Her- 
mitian operator belonging to different eigenvalues. Thus, 
we can use the Kroneker symbol to write 

rB 

U^Kr, Kb)Up(K r, n'b)rdr = 



J 2 M) 

JZJkB) 



- 1 



(Bll) 



Here, we sacrificed the beauty of symmetry and used the 
upper line of the Eq. (|B9|I : the same final answer is 
obtained from the lower line. 



b. Slightly less trivial part 

We have to perform now the limit B — > oo. The diffi- 
culty is that when B changes, so does also k, since it is 
subject to boundary condition (|B8I) . To circumvent this 
problem, the following trick is suggested. Let us choose 
some particular value of k, then boundary condition l|B8|) 
is satisfied by some discrete set of B values. Let us send 
B — > oo stepping over these specific values and thus keep- 
ing k fixed. Then, when B is already large enough, we 
can resort to the well known asymptotics 



Jfj,{x) ~ \j — cos 



2 

TTX 

2~ 



sin 



X ~ 2 lM 



Then, formula (|B11|) yields 

pB 

U^nr, nfyU^n'r, n'b)rdr 



ttkB 



jVKb) 



2 cos 2 [kB- § (n+ §)] 



(B12) 
(B13) 



On the other hand, we can also use asymptotics IB 12(1 for 
Jfj,(nB) and for Y^{kB) to simplify the boundary condi- 
tion (|B8|I : by some easy manipulations, we can re- write 
this boundary condition in the form 



1 



Y^nb) 



Jfj.(Kb) 



1 



(B14) 



" cos 2 [kB- § + 
Then, formula (|B13(I yields to the leading order in B: 

rB 

Ufi^nr, Kb)U^(n'r, n'b)rdr = 

= — [j 2 ( K b) + Y 2 ( K b)]A KK , (B15) 

Finally, we argue that at large B the Kronecker A 
should be replaced with Dirac S according to 



—6(k-k') 



(B16) 



To make this conclusion, we switch to the view point in 
which B can be arbitrary, while —k 2 and — k' 2 are the 
eigenvalues which depend on B. Then, when B — > oo, the 
eigenvalues come closer to one another, with the interval 
between neighboring k equal to w/B, as it is clear from 
the asymptotics of Bessel functions (IB12|I . Therefore, 
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any sum involving Kroneker A can be transformed into 
the integral 



dKf 



(B17) 



For completeness, we also mention the large k asymp- 
totics of Z, which turn out to be particularly nice: 



Zn(nr, Kb) ~ \ sin (n(b — r)) 

TTKT 



(B21) 



which means precisely l|B16ll . 

Taken together, equations (|B15|) and (|B16|) yield the 
answer 



U^Kr, k&)C/ m (/vV, n'b)rdr 



J?(Kb) + Y2(Kb) 



S(k-k') 
(B18) 

which is essentially formula IjBlfl . Formula (|B2|) follows 
automatically from (|B1(I and the the fact that functions 
Z^ form a complete set, which, in turn, follows from the 
very general spectral consideration. 



2. Asymptotics of Z 

Here, we first briefly describe the derivation of the 
small k asymptotics of Z^(kt, Kb), Eq. I|15|) . Knowing 
that El 



MO 



(£/2)" 

t(i + m) 



r(i + M ) 



T(l — ji) sin7r^ 



at £ 2 -c 1 + n, and using the identity 

r(l + ^)r(l - fi) sin 7T/z = Trfj, , (B20) 

and not resorting to any further approximations, we ar- 
rive at the Eq. (fT5|) . 



APPENDIX C: ALTERNATIVE 
REPRESENTATION OF THE GREEN'S 
FUNCTION 

Consider again the diffusion equation, which at t > t' 
reads 



dG a 2 



AG, 



dt 4 

subject to the initial and the boundary conditions 



(CI) 



G\ t=0 = 5(r-f) 
Let us denote 



This satisfies 



G\. 



, -> . (C2) 



o 

/•- Or dr 



+ OC 



e - P t / e ^ e G(t,fy)d9dt . (C3) 



1 d 2 8 



r 2 ^-G p ^(r : r') 



G P .A r ' r ') 



= ^5{r-r') 



(C4) 



The solution of this equation is the linear combinations 
of the Bessel functions and 1^ . Making it to satisfy 
the boundary conditions at r = 0, r — ► oo, and at r = r' 
(the latter dictated by the <5-function) , one arrives at 



G P Ar,r')=l 



(^) K, 



4 Kj^) 



2&VP 

r 



-I 



-I, 



(^) (^) 

(:-¥) (^) 



when r < 



when r > 



(C5) 



As expected, this is the symmetric function of r and r' ■ 
This formula was already obtained in [f| (equation (2.9) 
of that work). 

What we should do now is to invert the respective 
Laplace and Fourier transforms: 

G(r,r',6,t) = -\ r [°° I e-^ S e pt G P Ary)dpd», 

(C6) 



where C is the vertical contour in the plane of complex 
variable p which should be to the right of all singularities 
of G p ^ t . Knowing the explicit expression of G PiM , Eq. 
(IC5JI . we see that it has the singularity at p = 0. This 
singularity is due to both the branch point of y/p and 
the singular behavior of many Bessel functions at zero. 
Then, it is convenient to place the branch cut along the 
negative real axis in complex p-plane, and then to de- 
form the contour from C to C\, as shown in the Figure 
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QJ Then, because of the branch cut, on the lower side of 
the contour C\ we have p = e~ l7r \p\, while on the upper 
side we have p = e l7r \p\. Furthermore, instead of it 
is convenient to introduce the new variable, k, such that 



\p\ = K 2 a 2 /4. Then, integrals along the lower and along 
the upper sides of the contour C\ are each represented by 
integration from to oo over k. We can combine these 
two integrals together, and then simple algebra yields 



G(r,r',6,t) = |~ e-^e-M*/* { [I, (t«r) K, (i K b) - J M (t«6) ^ (ucr)] - 

- IF ( /~ ZKr ,? (-^0 K i* (-* Kb ) - h (-™ b ) K n } ndnd/i . (C7) 







branch cut 


: ® 






> 

Ci 






:c 

■ 



FIG. 4: Integration contours on the complex p-plane. Expla- 
nations are in the text. 



The expression in curly brackets here can be simplified 
using the following three relations: 

Ifj, (ikt) (inb) - 1^ (tub) (inr) = 
= 1^ (-iKr) (-mb) - 1^ (-1Kb) (-ikt) = 

= -^[J^(Kr)Y^(Kb)~J^Kb)Y^(Kr)} , (C8) 



and 
and 



(1Kb) (-1Kb) = (|) 2 [J 2 (Kb) + Y 2 (Kb)] , (C9) 



(ikt 1 ) (-1Kb) - Kfj_ (-ikt 1 ) (iKb) = 

2 

= ~ [Jm i Kr ') Y v ( Kb ) ~ j m (* b ) Y n M] ■ (CIO) 

Using these three results, we directly see that the formula 
(|C7|) gets transformed into J§J. This can, of course, be 
considered as another proof of normalization conditions 
for ^-functions, Eq. I|B1IB2|I . 



APPENDIX D: INTEGRATION OVER THE 
COORDINATE r' OF THE TRAJECTORY END 



consider the limit of large t in the sense that to 2 ^> r 2 ; 
that means, by the time t the walker should have trav- 
eled from its starting point typically much further than 
to the origin, O. Then, we note that although the inte- 
gration over r' runs to infinity, the integral is dominated 
by r' up to about a^/t, because of the truncation by the 
exponential factor e r ' a . Accordingly, the argument 
£ = 2rr' /a 2 t of 7 M in the Eq. JSJ is small, and we can 



use the expansion I M (£) ~ 
gration over r', this yields 



(£/2)7r(l + (i). Upon inte- 



W(0) 



cx e 



2 /a 2 t 



cos(fiO) 



r 1 + h 



r(i + M) 



2 \ 2 



r 



dp, . (Dl) 



Now, we have to remember that r 2 /a 2 t <C 1, which means 
that the latter integral is dominated by small fi, more 
specifically by p up to about 1/ |ln (r 2 /a 2 t)\. Replacing 
both T-functions with unity leads then to 



W(9) 



hi 



oc 



6» 2 



(D2) 



where the latter transformation is justified again because 
t is large. Thus, the resulting distribution is indeed inde- 
pendent of r, and it is nothing else but the Spitzer law, 
Eq. P. 



APPENDIX E: PROOF THAT W IS THE 
PROBABILITY 

Here, we check that W satisfies the normalization con- 
dition as the probability: 

°° / -i ' \ 



Most easily, integration in I|11J) can be addressed using 
Eq. Indeed, whatever is the value of r, we should 



which also means that identification of the points 9 and 
8 + 2-kti erases all the topological information. 
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n— — oo 



First, let us denote for brevity z = 2rr' /a 2 t, and then 
we write 



]T W (6 + 2im, z) = 

OO 

^2 cos ((6 + 2im)n) 

I — — OO 

oo 

V" cos ((0 + 27rn) /z) 



= 2e 



I fl (z)dn = 
I M (z)dfcE2) 



Here, integration is expanded over all /i, both positive 
and negative, the price being the absolute value of fi 
serving as an index of 7ui. The expression in the square 
brackets can be easily transformed using the identity 



J2 e a***t = J2 6(t-m) 



(E3) 



h— — oo 

Thus, we write 



^2 cos ((0 + 2nn) fi) = 

71 — — OO 
1 OO 

t r e ^(0+27rn) _|_ e i(i{-0-27ni) 

71 — — OO 
1 OO 

_ |" e *K0+27rrO _j_ e *M(-^+27rn) 

n— — oo 
1 oo 



-ip,0~\ 



COS 



cos({j,9) — m) 



(E4) 



We also use the two identities l34 



I, =e- w *' 2 J v {iz) 



(E5) 



(formula 8.406.3 in and 



Z2 COS 



JoO) + 2^Yj fe (z)cos(/c0) (E6) 



fe=i 



n— — oo 



(formula 8.511.4 in j3J|). This yields 

oo 

E W(6» + 27rn,z) = 

OO 

— z cos \ A im6 7- / \ 

m— — oo 

oo 

-;cos0 e™ e - l|m|7r/2 J| m |(iz) = 

m— — oo 

oo 

Jo(zz) + 2 cos(m0 — mir)i m J m (iz) 



.,(«) «*(«-*■) 



x e- 2 ™ 3 " = 1 



(E7) 



This completes the proof. 
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